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Abstract 

We develop the nonadiabatic polaron theory of superconductivity of M^C^q taking 
into account the polaron band narrowing and realistic electron-phonon and Coulomb 
interactions. We argue that the crossover from the BCS weak-coupling superconduc- 
tivity to the strong-coupling polaronic and bipolaronic superconductivity occurs at the 
BCS coupling constant A ~ 1 independent of the adiabatic ratio, and there is nothing 
"beyond" Migdal's theorem except small polarons for any realistic electron-phonon in- 
teraction. By the use of the polaronic-type function and the "exact" diagonalization in 
the truncated Hilbert space of vibrons ( "phonons" ) we calculate the ground state energy 
and the electron spectral density of the C^q molecule. This allows us to describe the pho- 
toemission spectrum of CgQ in a wide energy region and determine the electron-phonon 
interaction. The strongest coupling is found with the high-frequency pinch Ag2 mode 
and with the Frenkel exciton. We clarify the crucial role of high-frequency bosonic exci- 
tations in doped fullerenes which reduce the bare bandwidth and the Coulomb repulsion 
allowing the intermediate and low-frequency phonons to couple two small polarons in 
a Cooper pair. The Eliashberg-type equations are solved for low-frequency phonons. 
The value of the superconducting Tc, its pressure dependence and the isotope effect are 
found to be in a remarkable agreement with the available experimental data. 



PACS numbers 79.60.Bm, 63.20.-e, 71.38.+i 



1 Introduction 



Any phonon mediated superconductivity is a nonadiabatic phenomenon. Nonadiabaticity 
of electrons in metals measured by the ratio of the characteristic phonon frequency to the 
Fermi energy is normally small (< 10^^), so is their superconducting T^. High temperature 
superconductivity of doped fullerenes H seems to be due to large nonadiabaticity. The 
phonon frequencies are high, uj < 0.2eV while the bare Fermi energy is very low Ep < 0.2eV 
0. Tolmachev logarithm in the definition of the Coulomb pseudopotential fi* does not 
apply in this nonadiabatic case and the electro n-phonon coupling should be strong (A > 1) 
to overcome the Coulomb repulsion. The strong electron-phonon interaction implies the 
nonadiabatic polaron or bipolaron superconductivity as discussed in detail by Alexandrov and 
MottQ]. However, the final answer on the origin of high Tc depends not only on the adiabatic 
ratio uj/Ep but also on the coupling constants with different phonons . If a relatively weak 
coupling (A < 0.5) with low-frequency phonons dominates and the Coulomb repulsion is not 
very large, then the Migdal-Eliashberg theory is applied with the BCS ground state. On 
the other hand, if the coupling is strong and (or) high-frequency phonons are involved, the 
nonadiabatic polaron theory is more appropriate. 

In the past years several different calculations of the electron-phonon coupling constants 
have been reported for fullerenes[^ S 0' H El- Some of them yield the strongest coupling 
with the high frequency Hg modes with a moderate electron-phonon coupling, A < 0.5. On 
the other hand Picket et aZ0 predicted the strongest coupling with the high frequency Ag{2) 
mode and A ~ 3. The difference in calculated coupling constants is quite remarkable, and may 
result in a qualitatively different understanding of the nature of superconductivity of doped 
fullerenes. Therefore, the experimental determination of A and more extensive theoretical 
work are required. 

Recent photoemission spectroscopy (PES) of the molecule Cqq ||T0[ allows us to estimate 
the relative contribution of different phonons. While the variations analysis by Gunnarsson 
et al in the truncated space up to five phonons showed the strongest coupling with a low- 
frequency Hg mode, we found by applying the polaron theory of PES that the coupling with 
the high-frequency Ag2 mode dominates [|r^. The coupling constant with this mode appears 



at least by factor ~ 1.5 larger than with any other mode in a qualitative agreement with the 
tight-binding analysis by Picket et a/Q. This demonstrate that an estimate of the electron- 
phonon coupling constants from PES, the Raman and neutron scattering using the canonical 
theory of metals may be inconclusive, since doped fullerenes are not conventional metals to 
which Migdal's "theorem" can be applied. 

There is some confusion in the literature concerning the violation of Migdal's theorem. As 
was realized long time ago [|12| the Migdal theory of electron-phonon effects in the normal 
state and the Eliashberg theory which generalizes BCS theory to incorporate Migdal 
theory in the normal limit, provide a rigorous basis for understanding in the weak and 
intermediate coupling regime. By the use of the 1/A expansion technique ||15|| one can readily 
show that small polarons appear at A ~ 1 and thus Migdal's theorem fails. The breakdown 
of the Migdal-Eliashberg theory at A ~ 1 has nothing to do with the nonadibatic corrections 
due to the so-called "crossing" diagrams (vertex corrections). Small polarons are formed at 
A > 1 both in the initially adiabatic u <^ Ep or nonadiabatic u > Ep system irrespectively 
of the ratio co'/£'^||T6[. Their formation is the result of the broken translation symmetry as 
discussed in ref.|l^, Of course, in special cases Migdal's theorem will not hold even for 
small A < 1: first, if either characteristic phonon has small momentum q <^ qo {qo is the 
Debye momentum), and second, if the Fermi surface has a nested topology. In these textbook 
examples |l^ the "crossing" diagrams are no longer small. Several authors [^, ^ |22|, p3 



suggested a theoretical treatment of the vertex corrections with quite opposite conclusions 
about their role for Tc. While many authors generated the opinion that vertex corrections 



should, in general, be negative, Grimaldi et al p3[ obtain a strong enhancement of Tc from 
non-adiabatic terms. 

In this paper we develop the theory of the superconducting Tc which takes into account the 
nonadiabatic effects, strong coupling with some phonons, and the realistic Coulomb interac- 
tion. First, we discuss in more details the origin of the breakdown of Migdal's theorem due to 
the broken translation symmetry. We argue that while a strong enhancement of Tc due to ver- 
tex corrections [^] is an artifact of an unrealistic electron-phonon interaction, the polaronic 
band-narrowing provides such enhancement quite naturally. It means that the multiphonon 



dressing is important for high- Tc superconductivity as was predicted by one of us [E^. Then 



we analyze the PES data [IC] for Cqq using the polaronic displacement transformation for 



Ag2 mode and the truncated Hilbert space for other phonons. Our fit to the experimental 



PES is just as good as the variations approach |TI]] for the low binding energies and better 
than that for the high-energy region. We obtain the strongest coupling with the high-lying 
Ag{2) pinch mode and with the Frenkel-type excitons. As a result, we determine the bare cou- 
pling constants to formulate and solve the Eliashberg-type equations fully taking into account 
multiphonon dressing and nonadiabatic motion of carriers. These provide us with the value 
of Tc, its pressure dependence and with the isotope effect which are in excellent agreement 
with experiment. In conclusion we analyze a possibility for the bipolaron formation and 2e 
Bose-liquid ground state of doped fuUerenes. 

2 Broken translational symmetry and breakdown of 
Migdal's theorem. 

In doped semiconductors the carriers become small polarons or bipolarons at the intermediate 
value of the coupling constant A > 1 or A > 0.5, respectively [^. However, within the 
adiabatic Migdal description of electrons and phonons coupled by the linear electron-phonon 
interaction there is no instability at any value of A if the bare ionic plasmon mode is replaced 
by the acoustic phonon mode[0. The corrections to the normal state spectrum due to the 
coupling are adiabatically small (~ oj/Ep). In particular, the critical temperature of the 
BCS superconductor is adiabatically small, T^/ Ep < u/Ep <^ 1 for all relevant values of A. 
Therefore, the self-consistent Migdal-Eliashberg approach does not allow for the possibility to 
study the small- (bi)polaron formation in the intermediate and strong-coupling regime. This 
drawback is due to the basic assumption of the canonical theory that the Green's functions 
(GFs) are translationally invariant, thus G{r, r', t) = G{r—r', t). This assumption excludes the 
possibility of the local violation of the translational symmetry due to the lattice deformation 
followed by the self-trapping. To enable the electron to relax into the lowest polaron state, 
one can introduce an infinitesimal translationally noninvariant potential, which should be 
set equal zero only in the final solution for GFs As in the case of the off-diagonal 



superconducting order parameter, a small potential violating a translational symmetry drives 
the system into a new ground state at sufficiently large A. Setting it equal to zero in the 
solution to the equation of motion restore the translational symmetry but in a new polaronic 
band rather than in the electron one, which turns out to be an excited state. 

In the Holstein model, in which electrons interact with the local (molecular) phonons 
one can notice the polaronic instability of the Fermi liquid at A ~ 1 already with Migdal's 
diagrams for the electron self-energy. The interaction is described by the Hamiltonian 

He-ph = ^7(q)^(q)Ck+qCkC?q + H.C., (1) 



V2N 

where the coupling constant 7(q) and the phonon frequency ct;(q) are momentum independent. 
The electron self-energy S in the Migdal approximation contains two contributions, Em, 
Fig. la and E^, Fig. lb. Em — Ao; and therefore remains adiabatically small compared with the 
bandwidth 2D ~ A'"(0)~^ in the relevant region of the couphng (A < D/to), which guarantees 
the self- consistency of the approach. On the other hand for the molecular phonons E^ ~ D\n 
is not small and it turns out to be comparable or larger than the Fermi energy already at 
A ~ 1 for any filling of the band (n is the electron density per cell). As a rule this diagram, 
which is momentum and frequency independent, is included in the definition of the chemical 
potential Ep. While this is justified for a weak-coupling regime, E^ leads to an instability for 
a strong coupling. To show this let us consider a one-dimensional chain in the tight-binding 
approximation with the nearest-neighbor hopping integral D/2. The renormalized chemical 
potential is given by 

Ef = Dsin l ^^"^^ - 2D\n, (2) 

with A = 7^a;/2D. The system is stable if the derivative dEp/dn is positive, which yields the 
following region for the stabihty of the Migdal solution: 

For two and three-dimensional lattices the numerical coefficient is different, but the critical 
value of A remains of the order of unity. 

In doped semiconductors and metals phonons are screened and the diagram E^ containing 
only zero momentum phonons vanishes. Therefore, one can erroneously conclude that there 



is nothing to worry about as far as the apphcabihty of the Migdal approach. However, it is 
sufficient to violate the translation symmetry and then restore it to observe the instability of 



the bare band. The polaronic instability is essentially an adiabatic effect |jT6|, 18 1. Therefore, 
to see how the same self-energy leads to the polaronic collapse independent of the type 
of phonons we consider the extreme adiabatic limit of the classical deformation field 0(r) 
coupled with the electron field ^(r) as 

H = He + J drdr' [g{r - r')0(r'){V^^(r)V^(r) - n\r)} + h.c] + s^jV^WI^. (4) 

Here Hf. is the electron kinetic energy, s is the sound velocity, g{r) is the coupling constant 
with the Fourier component = i^q^^7(q) /V'2N, n^(r) is the periodic density of carriers 
respecting the translational symmetry, and = sq. Minimizing Eq.(4) with respect to the 
classical field 0*(r) we obtain 

s2vV(r) = J dv'g*{v - r'){n(r') - n°(r')}. (5) 

The solution is 

<^W = -Ei^e^">q-<}, (6) 

where riq is the Fourier component of the electron density n{Y) = (^/'^(r)^/'(r)). Substituting 
Eq.(6) into the Hamiltonian Eq.(4) we find that the adiabatic lattice deformation leads to the 
lowering of the electron energy by the value 

5Mr) = -2E^^e^n^q-<}- (7) 

q 

If the electron density is periodic, riq = the shift of the energy is zero. Otherwise, it is 
not. For example, one can consider a random statistically uncorrelated distribution with the 
ansamble average (riq^iq') = A^n^5q^_q/. In that example the chemical potential is shifted by 
the value 



^El7(q)lK 
q 



1/2 



2nXD, (8) 



which is practically the same as in Eq.(2). However, now 



Ai^ = i?p = ^El7(q)r^q (9) 



depends on the phonon spectrum, integrated over all Brillouin zone rather than on zero 
momentum phonons. The broken translational symmetry lowers the energy by the value 
2Ep ~ 2XD per particle. The corresponding increase of the deformation energy is Ep. There- 
fore, the system prefers to relax into the self-trapped state if Ep > D. The nonadiabatic 
vertex corrections neglected within the Migdal approach are also completely irrelevant for the 
small polaron formation and, therefore for the breakdown of this approach. The extension of 
the Migdal approximation to the strong-coupling limit is unacceptable because of the broken 
translation symmetry. 

3 Exact nonadiabatic solution and polaron band 
restoring translational symmetry 

In the small polaron regime, A > 1 the kinetic energy remains smaller than the interaction 
energy and a selfconsistent treatment of a many polaron system is possible with a '1/A' expan- 
sion technique in the nonadiabatic regime |T^. This possibility results from the fact, known 



for a long time, that there is an exact solution for a single electron without hopping term. 
Following Lang and Firsov |30| one can apply the canonical transformation to diagonalise the 



Hamiltonian in the site representation 

H=YlT{in- n)clcj +J2ujiq)id\d^ + 1/2) + a;(q)n4M,(q)rfq + H.c] + H„ (10) 

<ij> q q,i 

where T(m) is the bare hopping integral, i = (m, s), j = (n, s) {s stands for spin), Ui{q) = 
7(q)ea;p(2q • m)/-\/2iV, and He is the direct Coulomb repulsion. The diagonalisation is exact 
if T(m) = (or A = oo): 

H = e^He-^. ill) 



where 



S = Y.h,{ui{ci)d^-h.c.). (12) 



The electron operator transforms as 



Ci = Qexp ^^Mi(q)rfq - /i.c.j (13) 



and the phonon one as 

dci = dc^-Y^ niU*{q) (14) 

i 

From Eq.(14) it follows that the Lang-Firsov canonical transformation is the displacement 
transformation for the multi-polaron system shifting ions to the new equilibrium positions. 
In a more general sense it changes the boson vacuum. As a result: 

H = - f^^i,j)4cj - EpYl^i + ^Yl Vijclc]cjCi + t^q(rf]jrfq + 1/2) (15) 

i,j i i,j q 

where 

aij = T(m - n)6s.s'exp ^^[^^(q) - Uj{(i)]d^ - /i.c.j (16) 
is the new hopping integral depending on the phonon variables, 

Vij = V"c(m - n) - 2^t^q {ui{q)u*{q)^ (17) 

q 

is the polaron-polaron interaction comprising the direct Coulomb repulsion Vc and the at- 
traction via a nonretarded lattice deformation (second term in Eq.(17)). In the extreme 
strong-coupling limit A ^ oo one can neglect the hopping term of the transformed Hamil- 
tonian. The rest has analytically determined eigenstates and eigenvalues. The eigenstates 
|A^) = \ni,nq) are classified with the polaron 11^,3 and phonon riq occupation numbers and 
the energy levels are: 

E = (T(0) -Ep-^l)Yn^ + lYl ^-^^^ + E ^q(^q + V2) (18) 

i i,j q 

where rij = 0, 1 and riq = 0, 1,2,3, ....oo. The interaction term does not include the on-site 
interaction i = j for parallel spins because of the Pauli principle. The polaronic (Franck- 
Condon) level shift is expressed in a general form independent on the type of phonons in 
Eq.(9). 

Thus we conclude that the Hamiltonian Eq.(lO) in zero order of the hopping describes 
localized polarons and independent phonons which are vibrations of ions relative to new 
equilibrium positions depending on the polaron occupation numbers. The phonon frequencies 
remain unchanged in this limit |2^, |18|. The middle of the electronic band T(0) falls down by 
Ep as a result of a potential well produced by the lattice deformation due to the self-trapping. 



With the finite hopping term polarons tunnel in a narrow band because of the degeneracy of 
the zero order Hamiltonian relative the site position of a single polaron in a regular lattice. To 
show this self- consistently one can apply perturbation theory using 1/A as a small parameter. 
Because of the degeneracy terms of the first order in T(m) should be included in a zero order 
Hamiltonian Hq: 



Ho = ^{cr{m - n) - fi5ij)clcj + ^uj^{dl^dq + 1/2) 



where 



cr(m - n) = {{aij)) = T{m - n)6s,s'exp[-g^{m - n)] 



(19) 



(20) 



is the hopping integral averaged with the phonon equilibrium distribution. It is calculated by 
the use of the relation 

(21) 



For zero temperature one obtains 



a m — n 



\exp 



-EK*(q)-^Kq)]< 



exp 



^[Ui(q) -nj(q)]rfc 



where 



/(m-n) = -J2(\ui{q) 
^ q 



\u. 



(q)|^-2<(q)«,(q) . 



|0), (22) 



(23) 



The bracket is equalto unity for T = 0. The straightforward generalization for finite temper- 
atures [|^, yields 



9\^) = ^ E l7(q) rcot/i [1 - cos(q • m)]. (24) 

Earlier we have shown by the numerical diagonalization of a finite site Holstein model 
that the expression for the polaron bandwidth, Eq.(20) is valid for the nonadiabatic and inter- 
mediate regime t^; > D for all coupling strength. In the adiabatic case important corrections 
appear. However, the main exponential term remains almost the same in the strong coupling 
regime, A >> 1. As a result of this brief recourse to the multi-polaron theory, one can con- 
clude that the major effect of high-frequency phonons is the reduction of the bandwidth and 
of the Coulomb repulsion. While optical phonons cannot overscreen the repulsion they reduce 



its value significantly. Then the acoustical and and low-frequency molecular modes provide a 
net attraction. 



It was argued by one of us |2i| that the polaron band narrowing can provide a high value 
of Tc because of the enhancement of the density of states. This enhancement does not depend 
on the particular choice of the electron-phonon interaction contrary to the enhancement of 
Tc due to the vertex corrections [^. The latter can be obtained only if the electron-phonon 
scattering is dominated by small momentum transfer q qn- A simple estimate of the 
screening radius in doped fuUerenes and cuprates yields a value of the order or even less 
than the lattice constant. Therefore, a long-range electron-phonon interaction is ruled out by 



screening, which makes the positive effect of the vertex corrections quite unrealistic p9 

4 PES and bare electron-phonon interaction 

Because of the diversity in the theoretical results for electron-phonon interaction the coupling 
strength can be unambiguously determined only from the experiment. Allowing a small 
shift of the phonon frequencies due to screening the coupling with the intraball modes are 
expected to be the same in doped fuUerenes as in the C^^J^ molecule. The recent photoemission 
spectroscopy of a molecule Cqq [|101 allows us to estimate the relative contribution of different 
phonon modes and other bosonic excitations to the interaction. 

In this section we analyze the PES data [T^ using the exact polaronic diagonalization 



with respect to the Ag{2) mode and introducing the polaron— exciton coupling. We obtain the 
strongest coupling with the high-lying Ag{2) pinch mode and with a Frenkel-type exciton. As 
a result we provide a strong evidence for the nonadiabatic coupling with high-energy bosonic 
excitations in M^Cqq- 

The Hamiltonian at hand, describing three degenerate ti^ states coupled with phonons, is 
diagonalized with respect to the Ag2 coupling using the canonical Lang-Firsov displacement 
transformation as described above. The result is 

H = e'He-' = -E^^^ ^ ^t^^^ + ^ g'^oo, ^ ^^M^^^^ + ^ ^ (25) 

m=l u=l n,m=l f fJ.=l 

where Ep^^ = g'^ujAg2 is the polaron shift due to the Ag2 mode, 3x3 dimensionless matrix M 



is taken from ref . 1 30| 



+ VSQi V3Q2 

M = VSQi -VSQ^ + Q5 V3Q3 
\ V3Q2 VSQs 

and Uj^^fj^ are the phonon occupation numbers of eight five-degenerate Hg modes with the 

phonon operators Q'^ = dl ^ + d^^^. The interaction with Hg modes is responsible for the 



dynamic Jahn- Teller effect in Ceo- According to calculations singly ionised C^q is in 
the intermediate coupling regime, while the doubly and triple ionised molecules are in the 
strong coupling limit with respect to the coupling with Hg modes. Therefore a reasonable 
estimate of the ground state energy is obtained by taking into account only diagonal part of 
M. Nevertheless, to avoid any ambiguity we calculated the spectral function Ipoi{uj) of the 
Hamiltonian, Eq.(25) by the exact numerical diagonalisation in truncated Hilbert space (up 



to 4 phonons) for the Hg modes as described in ref. |^8|. A self-trapped exciton in neutral C( 



--60 



was observed and discussed by several authors [^, 0]. Because of the polaron-exciton 
coupling we add the same spectral function to the total spectral density shifted by the exciton 
energy ujex, and multiplied by the polaron-exciton coupling constant g"^^ as 

I{u;) = Ipol{uj) + gljpoli'^ + ^ex)- (26) 

This is an exact procedure if the interaction with excitons is linear as with phonons. Then 
we integrate I{uj) with the Gaussian instrumental resolution function of width ~ AlmeV [|T^ 
taking into account the damping of the exciton in the second (excitonic) contribution with 
the substitution uJex ~^ ^^ex + ilex- We thus can fit the PES in a wide energy region as shown 
in Fig.2 with g'^ being the fitting parameters (inset). The polaron-exciton coupling constant 
is found to be g'^^ = 0.5, the exciton energy uJex — 0.5eV in agreement with the luminescent 
data [^. The inverse exciton lifetime is estimated to be '^ex — 580cm~^. The coupling to the 
Ag{2) mode turns out most important in agreement with the tight-binding calculations 
and with the doping dependence of the phonon frequencies and the Raman linewidths. The 
Ag{2) mode shows a clear shift with doping from 1467cm~^ in undoped Cgo to 1446cm~^ in 
KsCqo and to 1431cm~^ in KqCqq. The linewidth increases by ~ 3cm~^ in the metallic state 
(K^Cqq) and comes back to the value of Cqq with further doping in an insulator KqCqq. At 



the same time high frequency Hg{7, 8) modes also show large broadening in the metallic state. 
However, NqxCqq, which does not exhibit metallic state with doping, shows the same strong 
line broadening and the bleaching of the Hg{7, 8) modes [^]. Consequently, the broadening of 
phonon lines does not provide information about the electron-phonon coupling. Moreover, as 
has been pointed by Gelfand |^ ^g(2) mode cannot decay into electron-hole pair, no mater 



how strong the electron-phonon coupling is. This decay is prohibited by the conservation of 
the energy because the estimated half-bandwidth is below the Ag{2) frequency. We believe 
that an increase of the linewidth of Hg modes in the metallic phase is associated with the 
crystal field splitting of the five-fold degenerate mode. If this splitting is not very large it may 
result in an increase of the Raman linewidth. Therefore, the large broadening of Hg modes 
does not contradict to our conclusion that their coupling with the electron is weak. 

In general, an estimate of the electron-phonon coupling constants from the Raman and 
neutron scattering using the canonical theory of metals may be inconclusive since doped 
fullerides are not conventional metals. In particular, the Migdal "theorem" can not be applied. 

5 Polaron theory of 

It follows from the PES analyses that the frequencies of essential bosonic excitations (phonons 
and excitons) strongly coupled with electrons {g ~ 0.6) are above or of the same order as 
the electron half-bandwidth in doped fullerenes. This fact as well as the observation of the 
phonon and exciton-sided bands in PES by itself favor the nonadiabatic small polaron theory 
0] rather than the adiabatic Migdal-Eliashberg approach to M^Cqo- On the other hand 
the total contribution of intermediate and low-frequency modes is not negligible, while their 
individual coupling is relatively weak {g < 0.3). This allows us to treat these modes within 
the Migdal-Eliashberg approximation by the use of the Eliashberg function 

1 ^ 

a^F{n) = - ^ X^ujJin - u,) (27) 

where = 1, ...8 is the sum over 8 Hg modes with 

20 

X. = --glu;,N{0), (28) 



and QuyUJu from the inset in Fig. 2. The additional coefficient 10 in the partial X^, appears due 
to the five-fold degeneracy and the nondiagonal terms in the matrix M as explained in ref. 
1^. At the same time the strong nonadiabatic interaction with the pinch mode is fully taken 
into account by the renormalization of the halfbandwidth D as 



r,2 



W = De-<^ (29) 

and by the reduction of the bare Coulomb repulsion as 

V, = V,-2g^ujA,,. (30) 

The excitonic contribution is taken into account via the high-frequency dielectric constant e. 
Because of the covalent nature of Ceo molecules we do not expect any significant dispersion of 
e, which is estimated as e ~ 4.4 [^. In a solid the exciton frequency appears to be rather 
close to the fundamental gap ~ 2eV. Therefore, we assume that the excitonic effect on the 
bandwidth is included in the LDA density of states (DOS) A^(0) = 1/2D. In the following we 
take into account the polaron band narrowing due to the coupling with the pinch mode only. 
A possibility remains that we overestimate this coupling at the expense of the coupling with 
the two high-frequency Hg modes as discussed by Gunnarsson et al If so, one should 



also include the contribution of these two modes in the polaron band narrowing and in the 
reduction of the Coulomb repulsion by a simple redifinition of the high-frequency coupling 
constant g. 

With g'^ujAg2 — O.OGel/ (inset Fig. 2) we obtain the lowering of the Coulomb repulsion 
by about Q.12eV . By the use of the Coulomb law, the dielectric constant e ~ 4.4 and the 
lattice constant a ~ lAA we estimate intracell bare repulsion as Vc — 0.25eV. This value 



is in agreement with the binding exciton energy in solid Ceo |]35|. As a result, the residual 



repulsion is estimated as O.lSel^. The polaron density of states determined as 

Np{0) = N{0)e<^' (31) 

appears to be about Np{0) ~ Estates / eV spin if the bare density A^(0) ^ Q. Estates / eV spin 
according to some LDA calculations for -f^sCeo 0. Because different LDA calculations 
yield rather different values of the bare A^(0) from 6.6 up to 12. 'bstates / eV spin the 



polaronic density can be in the range from 9states/eVspin up to about 17 states /eV spin. 
These values of the polaron density of states are compared with those measured with NMR 
{17 states /eV spin), thermopower {lAstates / eV spin) , and magnetic susceptibihty 
{lAstates / eV spin)^, and estimated from the heat capacity jump and ESR 
(~ lOstates / eV spin) for K^Cqq. We note that different from the ordinary metals the polaron 
heat capacity and the Pauli susceptibility are renormalized by the same amount in the pola- 
ronic systemj^. By taking 6 .Gstates / eV spin as the bare total DOS of three degenerate bands 
we estimate the bare intraband Coulomb pseudopotential as 

= In{0)Vc = 0.55, (32) 

and the renormalized one as 

1 

3 



/i* = -Np{0)V, = 0.4. (33) 



Now we are prepared to solve the Eliashberg equations for T,,: 

W 



poo 

Jo '((CJ„' -CJ„)2 + ^]2) ((^^, +^^)2 + ^2) 



2Q 



T 2 AM -1 
y ^ ; tan 



W 



X 



N{0) Jo 



dna^F{VL)- 



2n 



(34) 



First we show the normal DOS N{E) in Fig. 3 calculated with the analytical continuation 
of the normal state self-energy (~ Z — 1) to real frequencies. N{E) takes into account the 
finite bandwidth,the band narrowing effect and scattering by low-frequency phonons. these 
results in drastic changes of the renormalized DOS compared with the initial DOS. The latter 
is taken energy independent for simplicity. 



The result for Tc is 



(35) 



for K^Cqo, and 



T, = 35K, 



(36) 



for Rb^CQQ. These values are obtained with the Ehashberg function, Eq.(27) and with the 
bare (LDA) DOS A^(0) = 6.6states/eVspin and A^(0) = 1 .'b states /eV spin for K-^C^o and 
RbsC^o, respectively. 

Assuming the inverse square root dependence of all intraball frequencies on the carbon 
mass and the square root dependence of gl we calculate the isotope effect. Normally, in the 
polaronic superconductor the isotope exponent a = —din {Tc)/d\n (M) is negative because of 
the increase of the polaronic density of states with the increase of the ion mass M ( (7^ ~ v^M 
in Eq.(31)) [^. However, the usual positive contribution to a from low-frequency modes can 
change its sign, which is indeed the case here. By the use of the Eliashberg equations with 
the isotope mass dependent polaronic DOS, Eq.(31) we obtain 

a = 0.27 (37) 

for K^Cqo, and 

a = 0.31 (38) 

for RbsCQQ. These values are rather close to those reported in the literature 0]. 

Finally, we calculate the dependence of Tc on pressure P assuming that only hopping 
integrals depend on pressure. In that case the polaronic halfbandwidth W has the same 
pressure dependence as the bare one, which according to ref. is described by 

W{P) = W{0)exp (^^^ , (39) 

with the elastic modulus K ~ 2A.5GPa |^0[. Solving the Eliashberg equations with the 
pressure dependent bare and polaronic bandwidth we arrive with the curve of Fig.4, which 
agrees well with the available experimental data. 

6 Conclusion 

A simple estimate of the characteristic energy scale of phonons and electrons in M^Cqq clearly 
shows that the canonical Migdal- Eliashberg theory should be modified to include the polaronic 
nonadiabatic effects. As discussed in the beginning of the paper the polaron collapse of the 
band is the consequence of the broken translation symmetry and has nothing to do with 



the vertex corrections. The polaron band narrowing represents a universal origin of the 
high Tc value. We argue that the high-temperature superconductivity of doped fullerenes 
is an example of the nonadibatic polaronic superconductivity discussed by one of us earlier 
23 1 . In that regime the polaron halfbandwidth is below the characteristic phonon frequency 



but still above the effective attraction between polarons. Therefore, the canonical BCS- 
Migdal-Eliashberg theory can be applied with the finite polaronic bandwidth and reduced 
Coulomb repulsion. The high frequency bosonic excitations play an important role reducing 
the bandwidth and the repulsion. They can be integrated out with the canonical displacement 
transformation, while the interaction with the low-frequency modes is still described by the 
Eliashberg equations. As a result our theory gives the experimentally observed values of T^, 
the isotope and pressure effects in M^Cqq with the reahstic coupling constants based on the 
PES data, the Coulomb law and common sense. 

Our calculated polaron halfbandwidth is about W ^ 150meV, or less. Fig. 3. The total 
contribution of all vibration modes to the attraction is 2i?*°*"' ~ 330meV^ (insert in Fig. 2). As 
a result the attractive interaction between two small polarons is estimated as A = 2Ep — Vc — 
SOmeV. The condition > A, necessary for the Fermi-liquid approach seems to be satisfied. 
However, the system is very close to the bipolaronic instability A > W with the 2e Bose- 
liquid ground state. The final answer on the nature of the ground state remains with the 
experiment. The hallmarks of small bipolarons are those of small polarons, plus 

- superfluid phase transition similar to that of He^, 

- spin gap in the magnetic susceptibility, that is Xs ^ at T — if a singlet is the ground 
state, and the absence of the Hebel-Slichter peak in the NMR, 

- electrodynamics of the charged Coulomb Rose gas with the divergent Hc2(T), 

- double elementary charge 2e in the normal state. 

The first three features are clearly observed in doped cuprates M. The upward Hc2{T) 



curvature near Tc is measured also in M^C^o as discussed in Ref.pS 
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Figures 

Fig. 1 Self-energy in Migdal's approximation. 

Fig. 2 Polaron theory fit (full line) to the experimental PES (dashed line). Frequencies 
oj = ujy, coupling constants g = , and the contribution to the ground state energy E = 
for different modes are shown in the inset. For comparison we also show the coupling constants 
(g(Gun), inset) and the calculated variations PES (dotted line) of ref. |T^. 

Fig. 3 Polaron density of states in K^Cqq. The bare DOS is shown with the dashed line. 

Fig. 4 Pressure dependence of Tc for K^Cqq and Rh^Cao compared with the theory. 



